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A choice of first-order variables for the characteristic problem of the linearized Einstein equations 
is found which casts the system into manifestly well-posed form. The concept of well-posedness 
for characteristic problems invoked is that there exists an a priori estimate of the solution of the 
characteristic problem in terms of the data. The notion of manifest well-posedness consists of an 
algebraic criterion sufficient for the existence of the estimates, and is to characteristic problems as 
symmetric hyperbolicity is to Cauchy problems. Both notions have been made precise elsewere. 
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I. INTRODUCTION 

Well-posedness is a concept naturally associated with 
initial-value problems, whereby the norm of the solution 
at some time in the evolution is bounded by the norm 
of the initial data, independently of the data. The norm 
of the solution is defined by integrating on a time-slice a 
positive-definite quadratic expression in the fields. The 
norm of the data is the same positive-definite quadratic 
expression integrated on the initial slice. A well-posed 
initial-value problem turns out solutions that are stable 
with respect to small variations of the initial data, and is 
therefore naturally interpreted as a necessary condition 
for the existence of a stable numerical scheme (of the 
finite-difference kind) for the problem 1]. 

The well-posedness of the initial-value problem of the 
Einstein equations has received much attention since the 
pioneering work by Choquet-Bruhat '5| , and has turned 
out a rich spectrum of equivalent formulations of the Ein- 
stein equations with differing properties, such as number 
of variables, number and nature of the associated char- 
acteristic speeds and allowed gauge choices A strong 
motivation for such attention is the expectation that well- 
posed formulations might play a main role in extending 
the run time of numerical simulations. 

The understanding of the characteristic problem of 
the Einstein equations seems to have been evolving in 
the opposite direction. With no obvious issues of gauge 
and with a reduced number of fields, numerical simu- 
lations of spacetimes via the characteristic formulation 
have been performing at long running times with com- 
paratively high efficiency ^ [3 , and the attention seems 
to have centered around improved accuracy [^. But, be- 
cause little if anything is known about the stability of 
solutions of characteristic problems of partial differential 
equations, interest in the concept of a well-posed char- 
acteristic formulation of the Einstein equations has been 



rather limited 0- 

The main obstacle to the concept of well-posedness of 
a characteristic problem seems to be the peculiar split 
of the data set into a "normal" subset which is usually 
interpreted as initial data, and a "null" subset which has 
usually been thought to have no analog in the initial- 
value framework. The "null" subset has been interpreted 
sometimes as a fiow of information through the bound- 
ary of the region where the solution is being sought. So if 
one were to ask if there exists an a priori estimate for the 
solution of a characteristic problem, how should the an- 
swer even be framed? What is the information available a 
priori of the solution? The question is confusing because 
the null subset of data must be prescribed on a surface 
that is transverse to the characteristic initial slice and 
that contains the direction of evolution. One must then 
know such data variables at all evolution times, therefore 
it seems hardly appropriate to call it a priori informa- 
tion. But if the a priori information is confined to only 
the normal data, then it is impossible to obtain a truly a 
priori estimate because the solution will be affected by 
the "flow of information" . One must then include the 
null data somehow in the estimate, and one way to do so 
for the characteristic problem of the wave equation has 
been proposed by Balean and Bartnik 0i S 13 • 

But in essence, a priori information is any data that 
are known before the solution is known. There is really 
no issue of "evolution time" . If we allow ourselves to 
regard the null subset as a set of a priori information 
as good as the normal subset, then we have the same 
number of data as the initial value problem of the same 
equations, and the issue is simply to find a way to com- 
pare the size of the solution within the region of interest 
with the size of the data on both data surfaces. More 
precisely, if the n dependent variables of the solution of 
the problem are represented as a vector v, the m(< n) 
null data are represented by and the n — m normal 
data are represented by q'^, then the estimate would read 

\\v\\'<Ki\\wY + \\<lY) 
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where K would be a number that would not depend on 
(ui^iO^). The norms would have to be defined appropri- 
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ately in each space. Clearly an inequality of this kind 
would ensure that vanishing values of {w^,q'^) on their 
respective data surfaces would return a vanishing solu- 
tion V, and the size of small perturbations of the data 
(10*^,5°) on their respective data surfaces would control 



the variation of the solution. The term 



would be 



the integral of a quadratic expression in the null data on 
the corresponding data surface. Whether we refer to it 
as data or as a flow of energy is quite immaterial, but 
the parallel with initial value problems is maintained if 
we choose to call it data. 

With this viewpoint, a notion of well-posedness for 
generic linear homogeneous characteristic problems can 
be formulated in reference to the problem's almost- 
canonical form^3|- Even though the choice of coordi- 
nates is not relevant to the essence of a characteristic 
problem, let's assume that the equations are written in 
terms of coordinates adapted to the characteristic sur- 
faces, so we have (m, x, a;*) with i — 1 . . .r — 2 and with 
u labeling the characteristic surfaces. Then one can al- 
ways find a choice of dependent variables v = {w, q) with 
w — wi, . . . , Wm and q = qi, . . . , qn-m so that the system 
of equations takes the form 



N"q,u +N^q,^ +N'v,, +N\ 



(la) 
(lb) 



where the n — m square matrix A'^" is non-degenerate. We 
refer to this as the almost-canonical form of the char- 
acteristic problem (the canonical form has the identity 
matrix in the place of TV" and is achieved by the trans- 
formation q N"q, which always exists because TV" is 
non-degenerate). This system of equations requires the 
values of q to be prescribed at m = and the values of 
w to be prescribed at a; = in order for a unique so- 
lution to exist in the region to the future of u = and 
to the right of a; = 0, that is: a total of n functions of 
r — 1 variables ■ We will say that the problem is well 
posed if there exists a number k independent of {q^,w'^) 
and a positive-definite rt— dimensional symmetric matrix 
H such that 



\v\\'T<e^^{\\wY + \\lY) 



where 



\v\\t 



Ju+x=T 



(9°, 9°) 



M=0 



\w 



,0||2 



(w",w") 



(2) 

(3a) 
(3b) 
(3c) 



In each case, the notation ( , ) denotes the standard Eu- 
clidean scalar product of the appropriate dimension; e.g., 
{Wjw) = X^/Li'^/- The notion of well-posedness thus 
consists of the norm of the solution with respect to a 
spacelike surface u+x = T being bounded by the norm of 
the data, irrespectively of the data. The surface u+x = T 



intersects the data surfaces u = and x = 0, and the 
fixed value of T is arbitrary, but finite, and may be not 
expected to be large if the problem has non-constant co- 
efficients and non-vanishing undifferentiated terms. 

An equivalent representation of an almost-canonical 
characteristic problem is the form 



C''v,a+Dv = 0, 



(4) 



where a = u,x,i, the principal matrix C" has a block 
diagonal form diag(A^",0), the principal matrix has 
a block diagonal form diag(A^^,l), and the remaining 
principal matrices C" are completely arbitrary. We say 
that the characteristic problem is manifestly well posed 
if the following two conditions are satisfied: 

i) All the principal matrices C° are symmetric 

ii) The normal block of the principal-u matrix, A'^", is 
positive-definite and the normal block of the principal- a; 
matrix, A^^, is non-positive definite but such that A^" -I- 
A^^ is positive definite (that is, -A^" < A^^ < 0). 

If these conditions are satisfied, then one can take 
H = -\- and there exists k such that Q holds. 
The proof of this statement is developed elsewhere , 
and parallels an argument first used to determine a pri- 
ori estimates for the characteristic problem of the wave 
equation • The notion of manifest well-posedness of 
characteristic problems as defined here is a close analog 
of the notion of symmetric hyperbolicity for initial- value 
problems 0|. The conditions (i) and (ii) together con- 
stitute a sufficient algebraic criterion for well-posedness 
of characteristic problems. Clearly, the criterion is not 
a necessary condition (just as symmetric hyperbolicity is 
not necessary for well-posedness). 

In Section ^ we introduce the characteristic prob- 
lem of the Einstein equations in second-order form -the 
form in which the equations are most commonly used 
for numerical implementation- and linearize them around 
Schwarzschild spacetime. Subsequently, we find a set of 
first-order variables that cast the characteristic problem 
into manifestly well-posed form in the sense introduced in 
the previos paragraph. We conclude that the linearized 
characteristic problem of the Einstein equations admits 
a well-posed formulation. Section IIIII offers concluding 
remarks regarding the reach of the result and some re- 
maining open questions. 



II. THE CHARACTERISTIC PROBLEM IN 
GENERAL RELATIVITY 

A. Original problem 

The line element ds^ = guudu^ + gurdudr + g^adudx'^ + 
Qabdx'^dx^ in coordinates adapted to a slicing by null hy- 
persurfaces labeled by a coordinate u is 



ds' 



r^habWU'' ] du^ - 2e^'^dudr 



-ir'^habU^dudx" + r^hatdx^dx'' 



(5) 
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where x°' = {9, 4>) label null geodesies in the constant- 
u slice, and det{hab) = sin^ 9, so that r is the lumi- 
nosity distance and functions as a parameter along the 
geodesies. We define h"''' as the 2-inverse of hab, via 
h°'^hbc — 5°'c- Enough coordinate conditions have thus 
been imposed to leave only six metric functions to deter- 
mine by means of the Einstein equations. 

Sachs' approach ^13] to the Einstein equations adapted 
to a null slicing considers the components of the Einstein 
equations G^v = in a frame adapted to the null slicing. 
Such a frame has the direction ^" normal to the slicing 
as its main leg. With the current signature convention, 
(—,-!-, the second leg of the frame is another null 

vector rf such that t°'na = —1. The frame is completed 
by a complex null vector m" normalized by m^rha — 1 
that is perpendicular to both and Ua. The complex 
vector m" is equivalent to two real orthonormal space- 
like vectors e", tangent to the surface of fixed value of 
u, but the complex combination m" = (e" -I- )/V2 is 
usually preferred because it facilitates some calculations 
while leading to more compact expressions. Up to scal- 
ing, the normal direction is ^" = g^^Ujp — g"", and is 
tangent to the null geodesies of fixed value of x"" on each 
sHce, being thus proportional to d/dr. 

In the current notation, Sachs' approach consists in 
splitting the Einstein equations into three sets. The 
first set contains six main equations, which themselves 
split into the four hypersurface equations, 
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Gn^m^m'^ = 0, and the two main equations. 



G^i,mf^m'^ — 0. The second set has only one equation, 



the trivial equation, Gp^^l^n'^ = 0. The third set con- 
tains the three supplementary conditions, G,, 



G^j^ijn^n^ = 0. The six main equations are the equa- 
tions that determine the six unknown metric functions 
hab, P, V, ?7°. Sachs argues that by virtue of the Bianchi 
identities (V^C'^ = 0), the remaining equations can be 
thought of as constraints on the data for the unknown 
metric functions and can be ignored for our purposes. 

The main equations have been used for numerical sim- 
ulations for some timej^tarting with cases with high sym- 
metry [ll|lljll|llmi9] and moving to full general- 
ity by 1996 |20|"7see |2l| for citations and an up-to-date 
review of the numerical implementation of Sachs char- 
acteristic approach to the Einstein equations). To my 
knowledge, the hypersurface equations appeared explic- 
itly for the first time in taking the following form: 



— fi,r — -r;;h'^'^h^'^hab,rhcd.r 
16 



(6a) 



= {r^e-^^hacU^r).r - 2r^ (^Aa) + r^h-^D.had^r 

(6b) 

= 2Vr-e'^'^n + 2e'^^{D''Dcf3 + D''l3D^p) 
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re 



(6c) 



• ■ 2 
Here D is the covariant derivative associated with hab 



(i.e., such that Dchab = 0), and TZ is the curvature scalar 
of the 2-surfaces given by the intersection of the constant- 
u and constant-r hypersurfaces. The evolution equations 
as they appeared in are 

= m°-m''(^r{rhab,u).r - -^irVhab,r),r 

-2e'^{DaDbP + DaPDbP) + h,^aDb)ir^U'),r 



hachbdUWi + —habrDJJ-' 



+r^U''Dchab,r - r^D''U^bha)c,r + r^D^bWha^cr) 

(7) 

Since we are interested in a real (as opposed to complex) 
version of the evolution equations explicitly in terms of 
the metric variables, rather than implicitly through the 
vector m", we wish to "un-contract" the vectors m°'m!' . 
Given a generic tensor quantity Hab, the contraction 
m°'m^Hab — kills both its skew-symmetric part and 
its trace, so the "uncontracted" equivalent of {Tj) is sym- 
metric in a, b and traceless: 

= r{rhabM),r - ^{rVhab,r),r 

~2e^^{DaDbf3 + DaPDbP) + /lc(aA)(r't/'),r 
y-.4g 2/3 ^ j,'2, 



2 -hachbdKKr + jhab^rDcU"^ 
+r^U''Dchab,r - r'^D''U(bK)c.,r + r^D,^bU''ha)c,r 
+ -jhab i^^h'^''^ ,r{hcd,u — -^hcd^r 

+2e^f\D''Dc(3 + D'^pDcP) - D^{r^U'')^r 

+ ^—h,dKU'r] (8) 

In the following, we linearize the equations around a 
Schwarzschild background. In our current notation, 
Schwarzschild spacetime is given hy U"" = P = 0, V = 
r — 2m and hab =diag(l, sin^ 0) = qab- For spacetimes 
that are close to Schwarzschild we can define first-order 
variables hab and V via 



hab — Qab + hab 

V = r-2m + V. 



(9) 
(10) 



This implies that the covariant derivative Da is equal to 
the covariant derivative on the sphere Vq up to terms 
of first order. Keeping only first-order terms in hab, V , (3 
and in the equations one obtains their linearized ver- 
sion: 



r{rhab), 
(r — 2to) 



2m, 



{rhab),rr = {(rhab) ,r ~ hab) + 2V a^ bP 

2 r 



-'?afcV'=Vc/3-gc(aV,,)(r^[/'=),^ 



+ 29abV,(r^[/^),, 



(11a) 
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P.r = (lib) 

(lie) 



3 cos 6 ~ ~ 
2 sm 

In the proeess to obtain (|lld|l we used 
2 ~ 



(lid) 
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sm ( 



3 cos 0~ ~ 
+ 77^ — 7:1109,9 — hee 

together with the fact that 

q'^^'Kb = 

by virtue of det(/iab) =det(qah)- 



(12) 
(13) 



B. Well-posed first-order form of the linearized 
characteristic problem 

We define now a set of variables that casts the original 
system into first-order form, namely, such that no second- 
derivatives appear. Clearly there is an infinite number 
of ways to achive this purpose, one of which has been 
used previously 23] . However, the following choice is 
particularly convenient for our current purposes: 

Pab = {rhab).r, (14a) 

Qa = qac{r^U'),r-2P,a, (Mb) 
Ta = 13, a. (14c) 
Ja = V''hab + qac{r^U^),r. (Wd) 

It is clear that hab,r is encoded in Pab and that the 
set {f3a,U°^) is encoded in the set {Qa,Ta), since given 
{Qa,Ta) we can solve uniquely for {(3a, U°^). In terms of 
these variables, equations (|lla|l . (|llc(l and l|lld|l read 

(r — 2m) 

'^^ab,u ^ ^ab,r 

1 2m 

+V(b(3a) - T^Qab^'^Qc = {Pab - Kb) (15a) 

^ ' 2 r 

rQa.r+y'^Pab = 2rqacU^ ~ 4Ta + Ja {I5h) 

Vr + q'"'Tb,d - ^g'' Jb,d 

3cos9~ „ 
-Tr^K9,9-rS/cU'' 
2 smO 

+q'''(^q';fhab,d-q"'Kakb,c 



-q'"'n/hae,c) = 2/3 - hse - T^q^'KdQc 



+9"' (q'^Kd) ,cheb 
+g-(g'"^r^,),e/lae (15c) 



The quantities F^j^ that appear in the right-hand side of 
(|15c|) are the Christoffel symbols of the metric of the 2- 
sphere, qab, thus representing given functions of 0, but 
not of the dependent variables. As it stands, Eq. (|15c|) 
is sufficiently explicit for our purposes, as will be made 
clear shortly. Additionally, we have 



Ta,r = 
rJa,r = Qa-2Ta + 2rUa, 



(15d) 
(15e) 



as a consequence of (|llb|l and of (|llc(l . Furthermore, by 
definition, wc have the additional equations 

rhab,r = Pab - Kb (15f) 

r2[/° = Q° + 2r" - 2r[/". (15g) 

Equations (|llb|l, l(T5ajl . I|15b(l . (fT5c| . (|15d|l . ifT^ . 
(|15f|l and P5g| ) constitute a linear homogeneous 
characteristic system of equations for the variables 
{Kb,P,U'',V,Pab,Qa,Ta,Ja) in Canonical form. The 
two variables Pab are the normal variables of the prob- 
lem, whereas the remaining twelve variables are all null 
variables. The system has a unique solution for normal 

data P^"' = {qi,q^) = q^ given on the sficc u = 0, 

and null data (h^Hl P^^W^^^'W^^) ,Q^^\t^'\ /^^) ^ 
{wi, . . . , Wii) = given on the worldtube r = ro 11] . 

For our purposes, however, we point out the following 
fact, particular to the linearized regime of the charac- 
teristic equations. The variable V appears only in Eq. 
(|15c|) . The rest of the equations in the system yield the 
values of the remaining variables, which can be taken as 
a source for V in Eq. H15c|l . Thus the variable V can be 
thought of as decoupled from the system, in the linearized 
regime. In other words, the system can be solved for 
{Pab, Qa,Ta, Ja, hab, P, U"") and then the solutions can be 
used to integrate V from equation H15c|l . The subsystem 
that yields {Pab, Qa, Ta, Ja, hab, [3,11"^) is also a charac- 
teristic system in canonical form, with a two-dimensional 
normal space, and an eleven-dimensional null space. In 
the following, this is the point of view that we take in 
order to argue that the linearized system is well posed in 
the sense of Section 

The subsystem obtained by leaving out Eq. H15c|l has 
the form 



C"i; „ + Cv.r + Cv.a + Dv = 



(16) 



where v = {Pab,Qa,Ta, Ja,Kb, [3,U°-). The matrix C" 
has the form diag(l, 1, 0, . . . , 0), whereas the matrix C 
has the form diag(— (r — 2m,)/2r, 1, . . . , 1). Additionally, 
the matrices C° are symmetric, coupling only the vari- 
ables Qa and Pab- The system is linear and homogeneous 
and meets the criterion for well posedness outside of the 
event horizon, since the normal block of is negative 
definite for all r > 2m, and the normal block of C"-|-C"" is 
positive definite for all values of r. This implies that there 
exists a constant n independent of the data = (g°, w°) 
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such that ||w||^ < e'^^ + where the norms 

are defined as follows 

llwlll = [ w(C" + COw (17a) 

Ju+r=T 
Ju=0 

\\WT EE [ ^(^?)2 (17c) 

As argued in Section this in turn ensures that the so- 
lutions are stable with respect to small perturbations of 
the data. 

Notice that the surface u + r — T for fixed value of T 
is a spacelike surface with respect to the hyperbolic op- 
erator because C" + C is positive definite. Additionally, 
one can verify that the foliation by constant values of 
M -I- r is also spacelike with respect to the Schwarzschild 
metric, and thus it corresponds to some timelike coor- 
dinate for Schwarzschild space other than the standard 
Schwarzschild time. This is, of course, immaterial for 
the purpose of defining a positive-definite norm of the 
solution V. 



III. CONCLUDING REMARKS 

The main result is a formulation of the linearized 
characteristic problem of the Einstein equations that is 
a counterpart to a symmetric hyperbolic formulation of 
the 3-1-1 Einstein equations with respect to the stability 
of the solutions under small variations of the data. The 
equations are l|llb(l . (fT5ajl . I|15b|l . (|15d|l . ifT^ . ifTpl 
and | |15g| ), which are thought to be solved as a charac- 
teristic problem of partial differential equations, from 
normal data -P"j(?', a;°) on the surface u = and null data 
(QUu, x-),TSiu, J^{u, X-), h<i,{u, x'^), P°{u, Ugiu, 
on the worldtube r — rQ outside of the event horizon of 
the background Schwarzschild black hole. The solution 
to this problem is used as a known source for the 
remaining metric variable V which is obtained by inte- 
grating an ordinary differential equation in r, Eq. H15c|l . 
from data V^{u,x'^) given on the worldtube. The 
problem has 14 variables in all, a significant reduction 
with respect to the standard hyperbolic formulations 
of the Einstein equations which gcnerically require 30 
dependent variables. 

With respect to the reduction in the number of vari- 
ables, it is perhaps appropriate to point out that the 
reduction has not so much to do with the use of charac- 
teristic coordinates, but much to do with the fact that 
the original equations are of second-order. In this re- 
spect, consider the original second-order formulation of 
the Einstein equations in terms of one timelike and three 
spacelike coordinates: they are six second-order equa- 
tions for six variables where all second-order derivatives 
of all variables occur generically. The first-order repre- 
sentation of such a system generically requires that all 



first-derivatives of the six variables be promoted to new 
variables, that is: an addition of 6 x 4 new variables. 
Hence the generic number of variables of a first-order 
formulation of the Einstein equations: 30. Any one of 
such formulations is hyperbolic in the sense that it ad- 
mits characteristic surfaces: surfaces that are not ap- 
propriate as initial-data surfaces for the Cauchy prob- 
lem. One can write the characteristic problem of such 
a hyperbolic formulation, but this problem will have no 
less than 30 variables - the characteristic problem of any 
first-order hyperbolic system of equations has as many 
variables (and data) as the initial- value problem 0, [ll| . 

But now consider writing the second-order form of 
the Einstein equations in terms of a coordinate that la- 
bels characteristic surfaces of the second-order system. 
By construction, the second-order derivatives with re- 
spect to the characteristic coordinate do not occur in the 
equations (that's why the coordinate is characteristic). 
Therefore, the first-order representation of the charac- 
terisic problem of the Einstein equations must have at 
least six fewer variables than the generic first-order hy- 
perbolic formulation, that is, at most 24 variables. A 
simpler instance of the non-commutativity of the oper- 
ations of casting a system into characteristic form and 
of reducing to first-order is that of the two characteristic 
problems of the (scalar) wave equation. The problem of 
reducing to first-order and then going to a characteristic 
formulation can be found in 10], whereas the problem of 
first transforming to characteristic coordinates and then 
going to first-order appears, for instance, in |l2l |. 

The characteristic problem that we have examined 
here does not return a solution to the Einstein equations 
as it stands, because we have neglected four equations: 
the "trivial" equation and the "supplementary condi- 
tions" I which do not seem to appear in explicit form 
in the literature. Irrespectively of their explicit form, 
■at|:)ese equations are essential if one is to obtain a met- 
ric that solves the Einstein equations via the character- 
istic approach. Nevertheless, considering them jointly 
with the main equations for a global statement of well- 
posedness is outside of the scope of the current work. 

A number of issues remain open to further study. One 
question is whether estimates can be established for the 
derivatives of the solution of the characteristic problem of 
the Einstein equations. As stated in manifest well- 
posedness of a linear homogeneous characteristic problem 
in general does not guarantee the existence of a priori 
estimates for the derivatives (in stark contrast to the case 
of initial- value problems 24| ) . Another open issue is how 
the well-posedness of the characteristic problem relates to 
well-posedness of a hyperbolic formulation of the Einstein 
equations. 

Lastly, one may want to consider other first-order re- 
ductions of Sachs' characteristic problem in order to es- 
tablish estimates that are of greater relevance to the non- 
linear case, in particular the 18-variable first-order re- 
duction developed in On the other hand, whether 
or not the argument used in the linearized case is rele- 
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vant to the full non-linear case is not really an issue at 
this time, since a concept of well-posedness for non-linear 
characteristic problems is not available yet. 
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